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PURPOSE:  The  purpose  of  this  CHETN  is  to  examine  storm  surge  and  wave  reduction  benefits 
of  marsh  restoration,  as  well  as  the  impact  of  future  wetland  degradation  on  local  surge  and  wave 
conditions.  Storm  surge  simulations  of  representative  moderate  and  severe  hurricanes  were  per¬ 
formed  using  the  ADCIRC  storm  surge  model,  with  the  inclusion  of  radiation  stress  gradients 
from  the  STWAVE  nearshore  wave  model.  Coupled  model  simulations  were  made  for  wetland 
landscape  configurations  in  south  Louisiana  that  involved  hypothetical  restored  and  degraded 
wetland  features.  This  is  the  first  in  a  series  of  tech  notes  on  the  influence  of  marshes  on  storm 
surge  and  waves. 

INTRODUCTION:  Coastal  areas  are  vulnerable  to  devastating  storm  surge  and  waves,  a  vulner¬ 
ability  that  will  increase  with  the  ever-increasing  population  that  seeks  to  reside  along  the  coast 
and  its  accompanying  infrastructure.  The  threat  will  be  exacerbated  by  any  increase  in  the  fre¬ 
quency  and  severity  of  hurricanes  as  well  as  an  increase  in  relative  sea  level.  Topography,  land¬ 
scape  features,  and  vegetation  have  the  potential  to  reduce  storm  surge  elevations  and  absorb 
wave  energy.  Land  elevations  greater  than  the  storm  surge  elevation  act  as  a  physical  barrier  and 
create  bathymetric  resistance  for  the  surge  and  waves.  Landscape  features  such  as  wetlands  also 
have  the  potential  to  create  frictional  resistance  and  affect  storm  surge  and  wave  energy  even  if 
submerged.  The  purpose  of  this  study  is  to  apply  numerical  models  to  assess  the  potential  of 
coastal  landscape  features  for  reducing  storm  surge  and  waves  for  hurricanes  with  varying  inten¬ 
sity.  The  analysis  provides  valuable  information  on  trends  and  relative  performance  but  should 
not  be  taken  as  an  absolute  quantitative  assessment  of  surge  and  wave  reduction. 

MODELING  SYSTEM:  A  team  of  engineers  and  scientists  specializing  in  coastal  hydrodynam¬ 
ics,  meteorology,  statistics,  and  computer  science  collaborated  to  develop  an  integrated  modeling 
system  to  estimate  storm  surge  inundation.  A  schematic  diagram  of  the  system  is  shown  in  Fig¬ 
ure  1.  The  modeling  system  was  validated  for  Hurricanes  Katrina  and  Rita  (see  Bunya  et  al. 
2009)  and  was  shown  to  be  accurate  within  +/-0.5  m  of  the  measured  high  water  marks.  These 
storms  were  selected  due  to  the  unprecedented  quality  of  the  storm  system  data  and  resulting 
high  water  marks.  In  addition,  the  extent  of  inland  inundation  was  extreme,  allowing  for  a  unique 
opportunity  to  validate  the  effectiveness  of  numerical  models  that  include  variation  in  topogra¬ 
phy,  overland  resistance,  and  overland  wind  speeds. 

For  each  storm,  defined  by  a  track  and  time-varying  wind  field  parameters,  the  TC96  Planetary 
Boundary  Layer  (PBL)  model  (Thompson  and  Cardone  1996)  was  applied  to  construct  snapshots 
of  wind  and  atmospheric  pressure  fields  every  15  minutes  for  driving  surge  and  wave  models. 
ADCIRC  (Luettich  et  al.  1992,  Westerink  et  al.  1994,  Luettich  and  Westerink  2004)  was  then 
run  to  compute  the  pressure-  and  wind-driven  surge  component.  In  parallel  with  the  initial 
ADCIRC  runs,  the  large-domain,  discrete,  time-dependent  spectral  wave  model  WAM  (Komen 
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et  al.  1994)  was  run  to  calculate  directional  wave  speetra  that  serve  as  boundary  conditions  for 
the  local-domain,  near-coast  wave  model  STWAVE  (Smith  et  al.  2001,  Smith  2007).  Using  ini¬ 
tial  water  levels  from  ADCIRC,  winds  that  include  the  effects  of  sheltering  due  to  land  bounda¬ 
ries  and  reduction  due  to  land  roughness,  and  speetral  boundary  eonditions  from  WAM, 
STWAVE  was  run  to  produce  wave  fields  and  estimated  radiation  stress  fields.  The  radiation 
stress  fields  were  added  to  the  estimated  wind  stresses,  and  the  ADCIRC  model  was  run  again 
for  the  time  period  during  whieh  the  radiation  stresses  potentially  make  a  significant  contribution 
to  the  water  levels. 


Figure  1 .  Diagram  of  modeling  system  for  coastal  inundation  application. 


MODEL  SETUP:  The  modeling  system  was  applied  to  perform  a  sensitivity  analysis  to  assess 
the  impact  of  wetland  restoration  and  degradation  on  hurricane  surge  and  waves.  Simulations 
were  made  for  landseape  configurations  in  south  Louisiana  that  involved  a  base  condition  and 
hypothetical  restored  and  degraded  eonditions  for  Caernarvon  marsh.  The  base  configuration  was 
consistent  with  the  existing  post-Hurricane  Katrina  condition  in  south  Louisiana.  For  the 
degraded  condition,  the  marsh  area  outlined  in  Figure  2  was  lowered  to  0.6  m  below  North 
Atlantie  Vertical  Datum  1988  (NAVD88)  to  represent  an  estimate  of  the  long-term  effeets  of 
subsidence  and  degradation  and  charaeterized  as  open  water.  The  improved  marsh  was  raised  to 
approximately  +0.3  m  NAVD88  and  eharacterized  as  intermediate  marsh. 

Two  storms  were  simulated  for  each  configuration  giving  a  total  of  six  simulations.  The  first 
storm  (HURl)  was  a  large  hurricane  of  moderate  intensity  with  a  central  pressure  of  approxi¬ 
mately  960  mb,  a  radius  to  maximum  winds  of  39  km,  and  a  forward  speed  of  20  km/hr.  HURl 
parameters  are  similar  to  Hurricane  Hilda,  which  made  landfall  in  south  central  Louisiana  in 
1964.  The  second  storm  (HUR2)  was  a  large  hurricane  of  severe  intensity  with  a  eentral  pressure 
of  approximately  900  mb,  a  radius  to  maximum  winds  of  39  km,  and  a  forward  speed  of 
20  km/hr.  These  parameters  are  similar  to  Hurricane  Katrina  in  2005.  Both  storms  followed  a 
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track  oriented  primarily  south  to  north  through  the  Biloxi  marsh.  The  track  path  is  plotted  in  the 
base  condition  result  figures  (see  Figure  3  for  example). 
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Figure  2.  Study  area. 


The  landscape  degradation  and  restoration  were  represented  in  the  numerical  models  by  bathy¬ 
metric  and  frictional  resistance  changes.  The  base  condition  coastal  feature  landscape  land  cover 
type  was  taken  from  the  United  States  Geological  Survey  (USGS)  National  Land  Cover  Dataset 
(NLCD)  classification  raster  map  based  upon  Landsat  imagery  and  the  USGS  Landsat  Data  Gap 
study.  Each  NLCD/Gap  classification  had  an  associated  land  roughness  length  {zoiand)  and  Man¬ 
ning  n  value  (see  Westerink  et  al.  2008  and  Wamsley  et  al.  2009).  The  roughness  length  values 
were  obtained  from  the  Federal  Emergency  Management  Agency  (FEMA)  Hazards  U.S. 
(HAZUS)  study  and  are  defined  as  the  height  where  the  wind  speed  becomes  zero  using  a  loga¬ 
rithmic  wind  profile.  The  final  Manning  n  values  were  determined  through  the  hindcasts  of  hur¬ 
ricanes  Rita  and  Katrina.  All  Manning  n  values  were  within  the  expected  range  of  values  for 
their  land-use  categories  according  to  standard  hydraulic  texts.  The  roughness  length  and  Man¬ 
ning  n  values  used  for  the  Caernarvon  marsh  analysis  are  representative  of  an  intermediate  marsh 
with  Zoiand  =  0.11  m  and  Manning  «  =  0.04. 
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Figure  3.  Peak  surge  level  for  HUR1,  Base  condition. 

The  winds  input  to  the  ADCIRC  and  STWAVE  models  were  reduced  to  account  for  the  higher 
surface  roughness  through  a  directional  land  masking  procedure.  Briefly,  the  land  masking  pro¬ 
cedure  incorporates  the  distance  for  the  wind  boundary  layer  to  adjust  to  the  surface  roughness. 
To  represent  this  phenomenon  in  the  models,  wind  reduction  factors  were  computed  based  on  the 
weighted  average  of  roughness  coefficients  (zoiand)  within  10  km  in  the  upwind  direction 
employing  an  inverse  distance  weighting  function  (see  Westerink  et  al.  2008  and  Wamsley  et  al. 
2009). 

The  speed  at  which  a  storm  surge  propagates  (and  the  resulting  water  elevation)  is  affected  by 
coastal  landscape  features  through  bottom  friction  and  form  drag.  Bottom  friction  is  generated  by 
fluid  shear  stresses  at  the  seabed  and  flow-drag  resistance  is  generated  by  fluid  stresses  on 
objects  extending  through  the  water  column.  Bottom  friction  occurs  in  relatively  shallow  areas 
and  bottom  friction  and  flow-drag  resistance  can  occur  in  vegetated  areas.  The  ADCIRC  and 
STWAVE  models  presently  only  account  for  bottom  friction,  and  therefore  the  effect  of  form 
drag  can  only  be  approximated  by  increasing  the  bottom  friction  coefficient.  The  ADCIRC  and 
STWAVE  models  apply  a  Manning’ s-type  bottom  friction  formulation. 
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RESULTS:  Storm  surge  and  wave  reduetion  benefits  of  marsh  restoration  and  degradation  on 
loeal  surge  and  wave  conditions  were  examined  by  comparing  peak  water  level  and  significant 
wave  height  maps  for  each  of  the  altered  Caernarvon  marsh  configurations  to  the  base  condition 
for  both  storms  HURl  and  HUR2.  The  analysis  focused  on  changes  in  water  level  and  waves  at 
the  hurricane  protection  system  along  the  east  and  south  facing  levees  in  the  Plaquemines  and 
St.  Bernard  Parishes  (see  Figure  2).  As  would  be  expected,  restoration  and  degradation  of  marsh 
resulted  in  decreases  (for  restoration  cases)  and  increases  (for  degradation  cases)  in  both  surge 
and  wave  magnitudes  at  the  hurricane  protection  system.  In  general,  the  wave  change  patterns 
are  consistent  with  the  water  level  changes.  The  waves  may  be  controlled  by  increased  frictional 
resistance  in  the  shallower  water,  non-local  depth-limited  breaking  (e.g.,  at  the  edge  of  the 
marsh),  or  be  depth-limited  at  the  marsh. 

Base  Condition.  Figures  3-6  are  the  peak  water  level  and  wave  height  maps  for  storms  HURl  and 
HUR2.  The  maximum  offshore  wave  heights  are  11m  and  16  m  for  HURl  and  HUR2,  respec¬ 
tively.  Maximum  surges  for  HURl  in  the  pocket  between  the  south  St.  Bernard  and  Plaquemines 
levees  (Caernarvon  pocket)  are  2  to  2.5  m  and  maximum  wave  heights  are  less  than  0.5  m. 
Surges  and  wave  heights  increase  for  HUR2  with  peak  water  levels  of  2.5  to  5  m  and  maximum 
wave  heights  are  greater  than  HUR  1,  but  still  generally  less  than  0.5  m. 


Figure  4.  Peak  surge  level  for  HUR2,  Base  condition. 
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Figure  5.  Maximum  wave  heights  for  HUR1,  Base  condition. 


Figure  6.  Maximum  wave  heights  for  HUR2,  Base  condition. 
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Caernarvon  Marsh.  Figures  7-8  plot  the  difference  between  the  Caernarvon  marsh  restoration  and 
base  condition  peak  water  levels  and  maximum  wave  heights  for  both  HURl  and  HUR2.  The 
impact  of  restoration  of  the  entire  Caernarvon  marsh  was  generally  limited  to  the  Caernarvon 
area  for  both  HURl  and  HUR2.  The  change  in  bathymetry  and  frictional  resistance  slowed  surge 
propagation  resulting  in  an  increase  in  water  level  seaward  of  the  changes  and  a  decrease  in 
water  level  landward.  Maximum  surge  elevation  was  reduced  less  than  0.5  m  along  the  south 
St.  Bernard  levee  and  between  0.5  and  1  m  along  the  Plaquemines  levees  for  both  storms. 
Maximum  wave  heights  were  reduced  by  less  than  0.5  m. 

Figures  9-10  plot  the  difference  between  the  Caernarvon  marsh  degradation  and  base  condition 
peak  water  levels  and  maximum  wave  heights  for  both  HURl  and  HUR2.  The  marsh  degrada¬ 
tion  allowed  surge  to  propagate  more  rapidly,  resulting  in  a  rise  in  peak  water  level  relative  to  the 
base  condition  for  both  storms.  For  HURl,  peak  water  levels  increased  0.5  to  1.5  m  for  HURl 
and  by  as  much  as  about  2  m  for  HUR2.  Maximum  wave  heights  along  the  levees  in  the  Caer¬ 
narvon  pocket  area  increased  as  much  as  about  1  m  for  both  HURl  and  HUR2. 
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Figure  7.  Difference  in  peak  surge  level  between  the  restored  and  base  conditions. 
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Figure  8.  Difference  in  maximum  wave  height  between  the  restored  and  base  conditions. 


Figure  9.  Difference  in  peak  surge  level  between  the  degraded  and  base  conditions. 
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Figure  10.  Difference  in  maximum  wave  height  between  the  restored  and  base  conditions. 


SUMMARY:  The  purpose  of  this  analysis  was  to  apply  numerical  models  to  assess  the  potential 
of  marsh  for  reducing  storm  surge  and  waves  for  a  moderate  and  severe  hurricane.  The  impact  of 
landscape  features  such  as  wetlands  on  surge  propagation  is  a  relatively  new  application  for 
surge  and  wave  models  and  an  area  of  active  research  that  suffers  from  a  lack  of  quality  data. 
The  analysis  provides  valuable  information  on  trends  and  relative  performance  but  should  not  be 
taken  as  a  definitive  quantitative  assessment  of  surge  and  wave  reduction.  The  set  of  hypothetical 
hurricanes  considered  here  is  quite  small;  and  there  are  many  other  possible  combinations  of 
hurricane  track,  intensity,  size,  and  forward  speed.  Different  combinations  of  hurricane  parame¬ 
ters  are  expected  to  change  the  effects  of  landscape  restoration/degradation  on  storm  surge  and 
waves  in  different  ways.  Results  indicate  that  coastal  marsh  does  have  the  potential  to  reduce 
surge  and  wave  magnitudes.  Restoration  and  degradation  of  marsh  resulted  in  a  decrease  (for 
restoration  cases)  and  an  increase  (for  degradation  cases)  in  surge  elevations  and  wave  heights 
for  hurricanes  of  both  moderate  and  severe  intensity.  The  magnitude  of  change  was  generally 
greatest  for  the  more  intense  storm. 

ADDITIONAL  INFORMATION:  Questions  about  this  CHETN  can  be  addressed  to  Mr.  Ty  V. 
Wamsley  (601-634-2099,  email:  Ty.Wamsley@us.army.mil).  This  Technical  Note  should  be  ref¬ 
erenced  as  follows: 

Wamsley,  T.  V.,  M.  A.  Cialone,  J.  Westerink,  and  J.  M.  Smith.  2009.  Numerical 
modeling  system  to  simulate  influence  of  marsh  restoration  and  degradation  on 
storm  surge  and  waves.  Coastal  and  Hydraulics  Engineering  Technical  Note 
ERDC/CHL  CHETN-I-77,  U.S.  Army  Engineer  Research  and  Development 
Center,  Vicksburg,  MS.  http://chl.wes.army.mil/library/publications/chetn/ 
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